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Abstract 

We analyze the method for calculation of a coefficient of jet azimuthal anisotropy with- 
out reconstruction of the nuclear reaction plane considering the higher order correlators 
between the azimuthal position of jet axis and the angles of particles not incorporated 
in the jet. The reliability of this technique in the real physical situation under LHC 
conditions is illustrated. 
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1 Introduction 



Nowadays the strong interest is springing up to the investigations and measurements of az- 
imuthal correlations in ultrarelativistic heavy ion collisions (see, for instance, [1] and references 
therein). One of the main reasons is that the rescattering and energy loss of hard partons in the 
azimuthally non-symmetric volume of dense quark-gluon matter, created initially in the nuclear 
overlap zone in collisions with non-zero impact parameter, can result in the visible azimuthal 
anisotropy of high-pr hadrons at RHIC [1, 2, 3, 4]. 

Recent anisotropic flow data at RHIC [5, 6, 7] can be described well by hydrodynamical 
models for semi-central collisions and pt up to ~ 2 GeV/c (the elliptic flow coefficient v 2 
appears to be monotonously growing with increasing p? [8] in this case), while the majority 
of microscopical Monte-Carlo models underestimate the flow effects (see however [9]). The 
saturation and gradual decrease of v 2 at relatively large transverse momentum {p T > 2 GeV/c), 
predicted as a signature of strong partonic energy loss in a dense QCD plasma, seem now to be 
supported by the preliminary data extending up to p T — 10 GeV/c at RHIC. The interpolation 
between the low-pr relativistic hydrodynamics region and the high-p T pQ CD-computable region 
was evaluated in [4] . 

The initial gluon densities in Pb— Pb reactions at y/s NN = 5.5 TeV at the Large Hadron 
Collider (LHC) are expected to be significantly higher than at RHIC, implying stronger partonic 
energy loss. Moreover, since at LHC energies the inclusive cross section for hard jet production 
at Et ~ 100 GeV is large enough to study the impact parameter dependence of such pro- 
cesses [10], one can hope to observe the azimuthal anisotropy for hadronic jet itself [11, 12]. 
In particular, CMS experiment at LHC [13] will be able to provide both the jet reconstruction 
and adequate measurement of impact parameter of nuclear collision using calorimetric infor- 
mation [14]. In the case of jets, the methodical advantage of azimuthal observables is that one 
needs to reconstruct only azimuthal position of the jet without measuring its total energy. It 
can be done more easily and with high accuracy, while the reconstruction of jet energy is a 
more ambiguous problem [14]. However the measurement of jet production as a function of 
azimuthal angle requires event-by-event determination of the nuclear event plane based on the 
anisotropic flow analysis. 

The methods for elliptic flow analysis can be generally divided in two categories: two-particle 
methods suggested and summarized in works [15, 16, 17] and multi-particle methods [18, 19]. 
In two-particle methods the contribution of non-flow (non-geometric) correlations to the deter- 
mination of the elliptic flow coefficient v 2 is of the order of 1/^/Nq, where N is the measured 
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multiplicity. In multi-particle methods this contribution goes down typically as 1/N ' , i.e., 
smaller by a factor of the order of N^ 4 . Thus experimental techniques based on higher or- 
der cumulant analysis should be able in many cases to allow access to the smaller values of 
azimuthal particle anisotropy in comparison with two-particle methods, due to automatic elim- 
ination of the major non-flow many-particle correlations and the systematic errors originating 
from azimuthal asymmetry of the detector acceptance. This kind of analysis for particle flow 
has been already done by the STAR Collaboration at RHIC [20]. 

In our previous Letter [21] we proposed the method for measurement of jet azimuthal 
anisotropy coefficients without direct reconstruction of the event plane, and illustrated its reli- 
ability in a real experimental situation. This technique is based on the calculation of correlations 
between the azimuthal position of the jet axis and the angles of particles not incorporated in the 
jet, the azimuthal distribution of jets being described by the elliptic form. To improve our ap- 
proach, in the given paper we extend our analysis [21], considering the cumulant expansion [18] 
of multi-particle azimuthal correlations. 

2 Correlators versus the jet elliptic anisotropy 

Let us remind some features of our previous investigation [21]. We start from the essence of 
techniques [15, 16, 17] for measuring azimuthal elliptic anisotropy of particle distribution, which 
can be written in the form 



Knowing the nuclear reaction plane angle ipn allows one to determine the coefficient v 2 of 
azimuthal anisotropy of particle flow as an average (over particles) cosine of 2ip: 



However in the case when there are no other correlations of particles except those due to flow 
(or such other correlations can be neglected), the coefficient of azimuthal anisotropy can be 
determined using two-particle azimuthal correlator without the event plane angle ipR. 
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Let us consider now the event with high-p r jet (dijet) production, the distribution of jets over 
azimuthal angle relatively to the reaction plane being described well by the elliptic form [11], 

= -±- [1 + 2vf cos 2(y? - V*)] , Nf = J d V —j— , (4) 

where the coefficient of jet azimuthal anisotropy v 2 et is determined as an average over all events 
cosine of 2ip, 

1 } , _ , , dNi et 

dip 

One can calculate the correlator between the azimuthal position of jet axis ipjet 1 and the angles 
of particles, which are not incorporated in the jet(s). The value of this correlator is related to 
the elliptic coefficients v-i and v 2 et as 

1 } J dN jet dN 

(< cos2(v^ - <p) >) event = —— J dip jet J dp> cos2(^ et -<p) 

U —TV —7T 

1 J dN jet ■ 

= j^t J d( Pjet cos2(ip jet - ip R ) -j—— v 2 = v 3 2 v 2 . (6) 

Using Eq. (3) and intermediate result in Eq. (6) (after averaging over particles cos 2(ipj et — ip) 
reduces to v 2 cos2(<pj et — ipn)) we derive the formula for computing absolute value of the coef- 
ficient of jet azimuthal anisotropy (without reconstruction of sign of v 2 et ): 

v jet = I <cos2(v9 jef -y9) > \ ^ 
\\/< cos 2((pi — tp 2 ) > I 

V event 

This formula does not require the direct determination of reaction plane angle tpR. The brackets 
( ) represent the averaging over particles (not incorporated in the jet) in a given event, while 
the brackets ( ) event are the averaging over events. 

The formula (7) can be generalized by introducing as weights the particle momenta, 

jet I < Pr(y?)cos2(y? ief -y?) > \ , , 

v 2(p) ~ \ — / = / • \°) 

\y< Pti(<Pi) Pt2(<P2) cos2(v9i - v? 2 ) >/ event 
In this case the brackets ( ) denote the averaging over angles and transverse momenta of 
particles. The other modification of (8), 

vj« = ( , < ^)cos2( ( p jet - ( p)> \ , (9) 

\yj< ^l(^l) £ 2 (y?2)cOs2(v?l-V? 2 ) >/ cvcnt 

(Ei((pi) being energy deposit in a calorimetric segment i of position ipi) allows one using calori- 
metric measurements (9) for the determination of jet azimuthal anisotropy. 

lr The other possibility is to fix the azimuthal position of a leading particle in the jet. In this case calculating 
azimuthal correlations can provide the information about azimuthal anisotropy of high-pr particle spectrum. 
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3 Higher order correlators 



The main advantage of the higher order cumulant analysis is in the fact that, as argued in 
Ref. [18], if flow is larger than non-flow correlations, the contribution of the latter to v 2 extracted 
from higher order correlators is suppressed 2 by powers of particle multiplicity iVo in an event. 
Thus, for example, the fourth order cumulant for elliptic particle flow is defined as [18] 

c 2 [4] = < cos2(y?i + ip 2 - f 3 - f±) > 

— < COS 2 ((/?! — (/9 3 ) > < COs2((/2 2 — fi) > — < cos 2 (if i — if 4) > < cos2((/2 2 — fs) > ,(10) 

and in the case of existing only correlations with the reaction plane (i.e. the factorization of 
multi-particle distributions is held as in Eq. (3)) is equal to 

c 2 [4] = -v\. (11) 

If now one defines the coefficient v 2 of azimuthal anisotropy through two-particle correlator 

v 2 = \J< cos2(y? 1 - f 2 ) >, (12) 

then the contribution of non-flow correlations, as argued in [18], is of order 1/y/No. While their 
contribution to v 2 , extracted from the fourth order correlator 

v 2 = (-c 2 [A]fl\ (13) 

scales as l/N^ 4 , i.e. it is suppressed by an extra factor of Corresponding data analysis 

based on Eq. (10) and result (11) has been already carried out at STAR [20]. 

Now using the derivation of Eq. (7) and the result (11) it is straightforward to obtain the 
formula for calculation (measurement) of coefficient of jet azimuthal anisotropy through the 
higher order correlator, which is less sensitive to non-flow correlations: 

= ( (_ C2 [4])3/J ~ < cos2 (Vjet + fi-f2- fz) > 

+ < COS 2(f jet - f 2 ) > < COS 2 (fl~ fs) > 



+ < cos2(f jet - f 3 ) > < cos2(^! - <p 2 ) >}) . (14) 

/ event 

Here we stress once more that in the case of existing only correlations with the reaction plane, 
Eq. (14) together with Eq. (7) transforms into identity. The formula just derived can be, as in 
2 This can be essential under data analysis with not high enough multiplicity of particles in an event. 
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Sect. 2, generalized for calorimetric measurements of energy flows: 



= ( (-c 2( ^[4]) 3 / 4 ^ ~ < ^1(^1)^2(^2)^3(^3) cos2(y? jet + (fi! - (p 2 - <p 3 ) > 
+ < E 2 ((p 2 ) cos2(ip jet - ip 2 ) > < E 1 ((p 1 )E 3 ((p 3 ) cos2(v?i - <p 3 ) > 



where 



+ < E 3 (ip 3 ) cos 2(ip jet - ip 3 ) > < E 1 (ip 1 )E 2 (ip 2 )cos2(ip 1 - ip 2 ) >} ) ,(15) 

/ event 

C2(E)[4] = < Ei ((pi)E 2 ((p 2 )E 3 (ip 3 )E A (ip A ) cos2(v?i + ip 2 - ip 3 - tp A ) > 

- < E l (<p l )E 3 (<p 3 )cos2(<p l - <p 3 )> < ^2(^2)^4(^4) cos2(y? 2 - (Pa) > 

- < Ei(v9i)E 4 (v9 4 ) cos2(y?i - <f 4 ) > < E 2 (ip 2 )E 3 (ip 3 ) cos 2(ip 2 - ip 3 ) > . (16) 

In the case when the azimuthal position of jet axis correlates not only with the reaction 
plane, one can try to improve this technique using the multiple correlators of another form: 
averaging over not all events but selecting some their sub-events. For instance, one can consider 
sub-events 1 and 2, when jets are produced with the rapidity y > and y < 0. Then calculating 
correlator 



C2 6 *[4] = ( — 1 [ < cos2(y? ie t - if + (f) jet - 0) > 

cos2(y?i - ip 2 ) X cos 2(0i - <f> 2 ) > 

+ < cos2((p jet - Lp - <f) jet + 0) > 



- < cos 2(ip jet - ip) > < cos 2((f) jet -(/))>]) , (17) 

/ sub— event 1, 2 

we obtain that, if there are flow particle correlations only and the distribution of jets over 
azimuthal angles is described by the elliptic form (4) in every sub-event, it is equal to 

cf [4] = vt(y > 0) vi et (y < 0). (18) 

In Eq. (17) the angles ip are defined as the azimuthal angles of particles and jets in sub-event 
with y > 0, and — in sub-event with y < 0. Correspondingly the brackets < > represent 



the averaging over particles in sub-events 1, 2, while the brackets ( j are the 

\ / sub— event 1, 2 

averaging over these sub-events. The generalization of Eq. (17) for calorimetric measurements 
of energy flow is obvious (similar to Eqs. (9) and (15)). We do not write also this result specially 
as the examples of utilizing six- and other higher order correlators. 
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4 Non-flow correlations 



Here we discuss the influence of non-flow correlations 3 on the v 2 et determination. There are 
various sources of such correlations, among which minijet production [22], global momentum 
conservation [23, 24], resonance decays (in which the decay products are correlated), final state 
Coulomb, strong or quantum interactions [25]. We restrict our consideration to two-particle 
correlations only. It will be enough to illustrate the advantage in using higher order cumulants. 
In this case multi-particle distributions are not factorized again and instead of Eq. (3) we have 

dN dN dN cor . 



Iff dN 

c 2 [2] = < cos2(y?i - (p 2 ) > = — J difi J dip 2 cos2(y?i - ip 2 ) [— 



.2 1 + ^corhl 
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Nl J** J dlf2 d^dV 2 ' 
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Eq. (6) remains unchangeable and result (7) transforms into 



V>«[2] = ( < C0S2 ^-V ? )> \ = vi * 1 + A . (21) 

\J< cos 2(^-^2) >/ VI + ^M 2 

V x 7 event 

After some algebra the fourth order cumulant (10) reduces to 



M r: . eor ,4 \\ (22) 



V - = l diPl I" 11 " 2 l d(Pi COS2 ^+^-^3-^4) 



3 See also Appendix of work [18] and Ref. [22] 
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~k J** /^3sin2(^-c, 3 ) (23) 



= ( v cor) 2 ^ d ^° y r 2 is an even function of angular difference (if 1 — <p 2 )- 
The numerator in Eq. (14) is rewritten in the following form: 

Num = - < cos2(if jet + (fi - if 2 - fz) > 

+ < cos2(if jet - if 2 ) >< cos2(y?i - ip 3 ) > + < cos2(y? ie i - if 3 ) >< cos2(</?i - y? 2 ) > 
= - ^cos2(v^-<M 1+ 2V - /Vl + 



1 + 3A 

+2vl cos2(if jet -^ R ) 1 ^ Vcar {^ + SIN, (24) 
where the terms SIN are proportional to sin 2(<fj et — ipp) and vanishing after averaging over 

At the first glance it is hard to see the advantage in using higher order cumulants from 
Eqs. (21), (22) and (24). However, it is reasonable to suppose that the contribution of two- 
particle correlations to the normalization factor N 2 is small, A -C 1, while their "second Fourier 
harmonic" v~ or can be of the order of v\. Then all direct two-particle correlations v~ or are auto- 
matically canceled out 4 from Eqs. (22) and (24) in first order in A, but are survived in Eq. (21). 
The non-direct two-particle correlations i>+ r , v + r + are survived. But they are suppressed in the 
comparison with the direct correlations v ~ or (contributing to f^'PD due to the fact that J^^ 2 
is an even function of angular difference (ifi — if 2) only in most physically interesting cases [22]. 
Moreover, for the small-angle 5-like correlations (-^%- ~ cxp( ~ ( ^ 2)2/2(j2) , a -> 0) and for 

' & \d<pid(p2 V27TO- ' 

the large-angle oscillating ones ( ~ cos2(</?i — if 2)) the non-direct correlations 
are equal to zero. Then 



1 -C2Y±\yi* / event 



in this case, while 



<V\ ~ vt J , 1 , _ . (26) 



1 + v~ r /v 



Thus Eqs. (25) and (26) demonstrate the better accuracy of higher order cumulants explicitly. 
4 This is one of the main advantage of the cumulant expansion. 
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5 Discussion 



In order to illustrate the applicability of the method presented with regard for the real physical 
situation, we consider the following model (see Ref. [21] for details). 

The model 

The initial jet distributions in a nucleon-nucleon sub-collision at y/s = 5.5 TeV have been 
generated using PYTHIA_5.7 [26]. We simulated the rescattering and energy loss of jets in 
gluon-dominated plasma, created initially in the nuclear overlap zone in Pb— Pb collisions at 
different impact parameters. For details of this approach one can refer to Refs. [10, 11]. Es- 
sential for our consideration here is that in non-central collisions the azimuthal distribution of 
jets is approximated well by the elliptic form (4). In the model the coefficient of jet azimuthal 
anisotropy increases almost linearly with the growth of impact parameter b and becomes maxi- 
mum at b ~ 1.2Ra, where Ra is the nucleus radius. After that v 3 2 et drops rapidly with increasing 
b: this is the domain of impact parameter values, where the effect of decreasing energy loss 
due to reducing effective transverse size of the dense zone and initial energy density of the 
medium is crucial and not compensated anymore by stronger non-symmetry of the volume. 
The kinematical cuts on jet transverse energy and rapidity has been applied: E 3 ^ 1 > 100 GeV 
and \y jet \ < 1.5. After this the dijet event is superimposed on the Pb— Pb event containing 
anisotropic flow. 

Anisotropic flow was generated using the simple hydrodynamical Monte-Carlo code [27, 21] 
giving hadron (charged and neutral pion, kaon and proton) spectrum as a superposition of the 
thermal distribution and collective flow. To be definite, we fixed the following "freeze-out" 
parameters: the temperature Tf = 140 MeV, the collective longitudinal rapidity Y™ ax = 3 
and the collective transverse rapidity Y™ ax = 1. We set the Poisson multiplicity distribution 
and took into account the impact parameter dependence of multiplicity in a simple way, just 
suggesting that the mean multiplicity of particles is proportional to the nuclear overlap function. 
We also suggested [21] that the spatial ellipticity of the "freeze-out" region is directly related 
to the initial spatial ellipticity of the nuclear overlap zone. Such "scaling" allows one to avoid 
using additional parameters and, at the same time, results in the elliptic anisotropy of particle 
and energy flow due to the dependence of effective transverse size of the "freeze-out" region 
on the azimuthal angle of a "hadronic liquid" element. Obtained in such a way azimuthal 
distribution of particles is described well by the elliptic form (1) for the domain of reasonable 
impact parameter values. 

To be specific, we consider the geometry of CMS detector [13] at LHC. The central ("bar- 
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rel") part of the CMS calorimetric system covers the pseudo-rapidity region \r]\ < 1.5, the 
segmentation of electromagnetic and hadron calorimeters being Arj x A0 = 0.0174 x 0.0174 
and ArjX A0 = 0.0872x0.0872 respectively [13]. In order to reproduce roughly the experimental 
conditions (not including real detector effects, but just assuming calorimeter hermeticity), we 
applied Eqs. (9) and (15) to the energy deposition Ei(ipi) of generated particles, integrated over 
the rapidity in 72 segments (according to the number of segments in the hadron calorimeter: 
72 x 0.0872 = 2tt; % = 1, 72) covering full azimuth. 

Note that in the CMS heavy ion physics program, the modified sliding window-type jet 
finding algorithm has been developed to search for "jet-like" clusters above the average energy, 
and to subtract the background from the underlying event [14]. Strictly speaking, after jet 
extraction the background energy deposition in the calorimetric cells should be redefined and 
can appear to be not exactly equal to the initially generated one. However we neglect this effect 
here. In a real experimental situation, in order to avoid the influence of jet contribution on 
the particle flow, one can consider jets and particles incorporated in the energy flow analysis 
in different rapidity regions. 

Numerical results 

We have found [21] that the accuracy of v J 2 et determination from Eq. (9) is close to 100% 
for semi-central {b < Ra) collision and becomes significantly worse in very peripheral collision 
{b ~ 2Ra), wherein decreasing multiplicity and azimuthal anisotropy of the event results in 
relatively large fluctuations of energy deposition in each segment. 

In the given paper we test the efficiency of the higher order correlator (15) and have found, 
at first, that the results for v 3 ^ 1 obtained from Eqs. (9) and (15) are practically the same. This is 
explained by the fact that our simple Monte-Carlo event generator gives the elliptic anisotropy 
of energy flow, correlated with the reaction plane, but no correlations between energy deposition 
in the calorimeter segments. We can introduce such correlations at the calorimetric level "by 
hand", simply assuming that the probability of finding the energy in a segment i and the 
energy Ej in a segment j is proportional to EiEj(l + c^-), where the "correlation strength" 
may be, for example, proportional to 5^ (the short-range 5-like correlations) or cos2((/?j — ipj) 
(the long-range oscillating correlations). In this case we became convinced that the higher 
order cumulant (15) was almost independent of such correlations (as it was shown in Sect. 
4), while the result of calculation (9) changed, closely following the formula (26) corrected by 
autocorrelation terms which are non- vanishing in finite summation [18]. 

We have also found at the calorimetric level that, taking into account the effect of possible 
detector inefficiency (i.e. that the particles and jets are not detected in a "blind" azimuthal 



9 



sector of size a), the accuracy of v 2 e determination appears to be less than 50% at a > 30° in 
our model calculation without correlations and at b > Ra, whichever algorithm (9) or (15) we 
used. 

Fig.l is presented to illustrate the improvement due to the fourth order cumulant method 
in the determination of jet azimuthal anisotropy v 2 et depending on the ratio v~ or /v 2 , where v 2 
is the coefficient of elliptic azimuthal anisotropy of energy flow defined here as 

V2 = 9 F , F > ( 27 ) 

and E max (i) and E min ^ are the maximum and minimum energy deposit in a segment respectively 
(i = 1,...,72). The coefficient v~ or determines the "correlation strength" at the calorimetric 
level, Cij = 72v~ or 5ij for short-range correlations 5 (the similar result is obtained for long- 
range correlations with c^- = 2v~ or cos2(y9j — ipj)). The plots show the 6-dependence of the 
"theoretical" value of v 2 et (calculated including collisional and radiative energy loss when the 
reaction plane angle is known in each event), and the v 2 et determined by the methods (9) 
and (15) for the three ratios v~ or /v 2 = 0, 0.01, 0.1. We used two values of the input 
parameter, the number of charged particles per unit rapidity at y — in central Pb— Pb 
collisions: dN^/dy = 3000 (Fig. la) and 6000 (Fig. lb). 

One can see that improvement due to the fourth order cumulant method (result of (15) is 
independent of v~ or /v 2 and coincides with the result of (9) for v~ or /v 2 = 0) is pronounced for 
more peripheral collisions, smaller particle multiplicities and larger "correlation strengths". 



5 Here one should note that the majority of sources of non-flow correlations mentioned above is effective at 
small angles between particles. In our case these correlations can be partially smoothed out after summing 
particle energies over the azimuthal angles in a calorimeter segment of finite size (~ 5°). This can result in the 
smaller value of the ratio v~ or jv\ in comparison with the ratio v~ or /v% (and, as consequence, in a somewhat 
less improvement due to the higher order method at the calorimetric level (15) in comparison with the particle 
level (14)). We still have no any adequate Monte-Carlo generator for particle flow effects at LHC including the 
physical model for correlations. Thus we can not estimate the real value of v~ or jv\ (true benefit from the higher 
order method) and just treat it here as a phenomenological parameter. 
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6 Conclusions 



In the present paper we have analyzed the method for measurements of jet azimuthal anisotropy 
coefficients without reconstruction of the event plane considering the higher order correlators 
between the azimuthal position of the jet axis and the angles of particles not incorporated 
in the jet. The method is generalized by introducing as weights the particle momenta or the 
energy deposit in the calorimeter segments. In the latter case, we have illustrated its reliability 
in the real physical situation under LHC conditions. Introducing in the model correlations 
between energy deposits in the calorimeter segments does not practically change the accuracy 
of the method using fourth order cumulant calculations (15), while the result obtained with 
second order correlator (9) is dependent significantly of the "strength" of such correlations. The 
advantage of the higher order cumulant analysis is pronounced for more peripheral collisions 
and smaller particle multiplicities. 

To summarize, we believe that the present technique may be useful investigating azimuthal 
anisotropy of jets and high-p r particles in heavy ion collisions at RHIC and LHC. 

Acknowledgments. Discussions with A.I. Demianov, S.V. Petrushanko, V.V. Uzhinskii, 
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Figure 1: The impact parameter dependence of "theoretical" value of v 3 ? including collisional 
and radiative energy loss (solid curve), and v 3 2 et determined by the method (9) for v~ or jv\ = 
(dashed curve), 0.01 (dotted curve) and 0.1 (dash-dotted curve). The result obtained using the 
fourth order cumulant method (15) coincides with the dashed curve. dN ± /dy(y = 0, b = 0) = 
3000 (a) and 6000 (b). 
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